arXiv: 1509.06725vl [math.CV] 22 Sep 2015 


EQUIDISTRIBUTION AND /3 ENSEMBLES 


TOM CARROLL, JORDI MARZO, XAVIER MASSANEDA, 
AND JOAQUIM ORTEGA-CERDA 


Abstract. We find the precise rate at which the empirical measure 
associated to a /?-ensemble converges to its limiting measure. In our 
setting the /3-ensemble is a random point process on a compact com¬ 
plex manifolds distributed according to the /3 power of a determinant 
of sections in a positive line bundle. A particular case is the spherical 
ensemble of generalized random eigenvalues of pairs of matrices with 
independent identically distributed Gaussian entries. 


1. Background and setting 

Let {X,u}) be a n-dimensional compact complex manifold endowed with 
a smooth Hermitian metric uj. Let {L, (p) be a holomorphic line bundle with 
a postive Hermitian metric (p. This has to be understood as a collection 
of smooth functions (pi defined in trivializing neighborhoods Ui of the line 
bundle. If ei{x) is a frame in [/*, then ||ei(x)|p = Thus (pi must 

satisfy the compatibilty condition (pi — cpj = log \gij\, where gij are the 
transition functions. 

As usual we denote by the global holomorphic sections. If L is 

a line bundle over X and M is a line bundle over Y, we denote by L Kl M 
the line bundle over the product manifold X x Y defined as L Kl M = 
'k*x{L) (g) 7ry(M), where vrx : X xY — )■ A is the projection onto the first 
factor and vry : X xY —)■ T is the projection onto the second. The line 
bundle LM M carries a metric induced by that of L and M. 

Given a basis si,..., sw of H^{L, X) we define det(sj(xj)) as a section of 
over by the identities det(sj(xj)) = sgn((T) 0^^ 

We fix a probability measure on X, given by the normalized volume form 
uj"', that we denote by a. 

Definition 1. Let pt > 0. A (3-ensemhle is an N point random process on 
X whieh has joint distribution given by 

(1) ^^11 da{xi) (g) • • • (g) da{xN), 

Zn 

where Zx = Zx{l3) is ehosen so that this is a probability distribution in X^. 

Observe that the random point process is independent of the choice of 
basis Sj. 

A particularly interesting case is when j3 = 2, since then the process 
is determinantal. Let K denote the Bergman kernel of the Hilbert space 
H^{L,X) endowed with the norm ||s|| = ||s(a;)||^(iiT(x). Here || • ||^ 
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denotes the norm induced by the metric 4> (see Section 1.2 for more details). 
Then 

II det(sj(a:j))|p = || det{K{xi,Xj))\\. 

Another interesting situation occurs when /3 —>■ oo. In this case the prob¬ 
ability charges the maxima of the function || det(si(xj))||. A set of points 
{xj}j with cardinality dim.H^{L,X) and maximizing this determinant is 
known as a Fekete sequence. The distribution of these sequences has been 
studied in [LOCIO 


MOCIO] , and [BB WN ll| and we will draw some ideas 


from there to study general /3-ensembles. 

We consider now the situation where we replace L by a power L^, /c G N, 
and let k tend to infinity. We denote by Nk the dimension of 
It is well-known, by the Riemann-Roch theorem and the Kodaira vanishing 
theorem, that 


dim/^°(L^A) = + ~ fe", 

n! 

where ci(L) denotes the first Chern class of L. 

For each k we consider a collection of Nj- points chosen randomly accord¬ 
ing to the law (0. For each k the collection is picked independently of the 
previous ones. 

Given points x[^\ ..., x^^j^ chosen according to Q, consider its associ¬ 
ated empirical measure ■ For convenience we will drop the 

superindex (k) hereafter. We are interested in understanding the limiting 
distribution of the measures /ifc. 

The following result is well known; see [BBWNil . 


Theorem (Berman, Boucksom, Witt). Let be the empirical measure 
associated to a Fekete sequence for the bundle X). Then, as k ^ oo, 


hk 


V : = 


{iddcf)"^ 

fx(idd(f>y 


in the weak-* topology. 


The measure n is called the equilibrium measure. 

There is a counterpart of this result for empirical measures of general 
/3-ensembles (see [Berl4 , which gives an estimate for the large deviations of 
the empirical measure from the equilibrium measure). Our aim is to obtain 
a different quantitative version of the weak convergence of the empirical 
measure to the equilibrium measure, measured in terms of the Kantorovich- 
Wasserstein distance between mesaures. 

This sort of quantification has also been studied, with different tools, in 

[mmT^ 


the context of random matrix models, (see for instance 


Med 3 


and MM14 , where similar determinantal point processes arise) 


In fact some of the /3-ensembles we are considering admit random matrix 
models, at least when dimc(M) = 1. For instance, Krishnapur studied in 
|Kri09 the following point process: let A, B he k x k random matrices with 
i.i.d. complex Gaussian entries. He proved that the generalized eigenvalues 
associated with the pair {A,B), i.e. the eigenvalues of A~^B, have joint 
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probability density: 


( 2 ) 


1 

Zu 


nir 


1 




ni 


1=1' ■ ' ' i<j 

with respect to the Lebesgue measure in the plane. 

It was also observed in Kri09 that, using the stereographic projection 

TT : §2 ^ C 
Pj 

the joint density (§ (with respect to the product area measure in the product 
of spheres) is 

^ n ~ -^iiiR3- 

^ i<j 

Since this is invariant under rotations of the sphere, the point process is 
called the spherical ensemble. 

A point process with this law had been considered earlier - without a 
random matrix model - by Caillol CaiSl as the model of one-component 
plasma. 

One typical instance of the process is as in the picture. 


Ensemble with 500 points 


Ensemble with 1200 points 


Ensemble with 2000 points 



The spherical ensemble has received much attention. We mention a cou¬ 
ple of properties related to our results. In [Bor 11 , Bordenave proves the 
universality of the spectral distribution of the k x A:-matrix A~^B with re¬ 
spect to other i.i.d. random distribution of entries. As an outcome, he proves 
that the weak-* limit of the spectral measures /ifc = ^ , where Xi are 

the generalized eigenvalues is the normalized area measure in the sphere. 
This convergence is rather uniform: in |AZ15 Alishahi and Sadegh Zamani 
estimate the discrepancy of the empirical measure with respect to its limit 
and give precise estimates of the Newtonian and the logarithmic energies. 


1.1. The Kantorovich-Wasserstein distance. To measure the unifor¬ 
mity and speed of convergence of the empirical measures /ifc to the limiting 
measure v we use the Kantorovich-Wasserstein distance W. Given proba¬ 
bility measures fj. and zz, it is dehned as 


W{fi, v) = inf 
p 



d{x,y)dp{x,y), 


where d{x, y) is the distance associated to the metric oj and the infimum is 
taken over all admissible transport plans p, i.e., all probability measures in 
X X X with marginal measures p and v respectively. 
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In general, the Kantorovich-Wasserstein distance is defined on probabil¬ 
ity measures over a compact metric space X, and it metrizes the weak-* 
convergence of measures. 

It was observed in LOC10| that in the definition of W it is possible to 
enlarge the class of admissible transport plans to complex measures p that 
have marginals p and v respectively. We include the argument for the sake 
of completness. 

Let 

(3) IK(/i, z^) = inf // d{x,y)\dp{x,y)\, 

J JxxX 

where the inhmum is now taken over the set S of all complex measures p on 
X X X with marginals p{-,X) = p and p{X, •) = v. 

In order to see that W{p, i^) = W{p, u), we recall the dual formulation of 
W (see |Vil09 Formula (6.3)]): 


(4) 


W{p,v) = sup I J^fd{p : / e Lipi i(X)| , 


where Lip^ i(^) is the collection of all functions f on X satisfying |/(x) — 
f{y)\ < d{x,y). 

For any complex measure p with marginals p and i' and any / G Lipj^ ;^(X) 
we have 


I fd{p-v) 

X 

Hence 


f{y))dp{x,y) 


XxX 


< 


d{x,y)\dp{x,y)\. 


XxX 


Wip,!^) < inf / / d{x,y) \dp{x,y)\ = W{p,i^). 

P^^JJxxX 

The remaining inequality {W{p,i') < W{p,iy')) is trivial. 

A standard reference for basic facts on Kantorovich-Wasserstein distances 
is the book [Vil09 . 

1.2. Lagrange sections. We hx now a basis of sections si,... of H^{L^,X). 
Given any collection of points (xi,..., xn^) we define the Lagrange sections 
informally as: 


ijix) = 


Sl(xi) ••• Sl(x) ■■■ Sl(xjv,) 


«iVi.(3^l) ••• SN,.{x) ■■■ SN^.ixNu) 


Sl{x-l_) 


si{xj) ■■■ si{xnJ 


‘^Ni.ixi) ■■■ SM,.{xj) ■■■ sjVj, ) 

Clearly ij G H^{L^,X) and ij{xi) = 0 if i / j and ||I'j(a;j)|| = 1. 

More formally, we proceed as in [LOCIO] : if ej{x) is a frame in a neigh¬ 
borhood Uj of the point xj, then the sections Si{x) are represented on each 
Uj by scalar functions fij such that Si(x) = fij{x)ej{x). Similarly, the met¬ 
ric k(j) is represented on Uj by a smooth real-valued function k(l)j such that 

To construct the Lagrange sections we denote by A the matrix 


(e-|0.GH/^^.(^^.)) 
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and define 

, Nk 

lj{x) := ^ijSi{x)^ 

where Aij is the determinant of the snbmatrix obtained from A by removing 
the i-th row and j-th column. Clearly Ij G H^{L^,X), and it is not difficult 
to check that ||^j(xi)|| = 5ij, 1 < < N}.. 

Notice that if we denote by Pkixi ,..., xat^,) = || det Sj(xj)|| then 


(5) 


Pk V^l 5 • ■ • 5 5 • • • 5 ) 


and thus E(||£j ||^) < 1. In the case of the Fekete points (/3 
1 by definition. 


oo), supx ||^i(a:)|| 


2. Main result 


Theorem 1. Consider the empirical measure pk associated to the fd-ensemble 

be the equilibrium measure. Then 


given in 


Definition\l\ and let 




_ (iddcpY' 

S^(idd4>Y 


mv{pk,i^) = o{i/Vk). 


Remark. The rate of convergence cannot he improved. Let a be any nowhere 
vanishing smooth probability distribution on X. Let be any discrete set 
on X with cardinality ffEk ~ fc"" ~ iV^, and let pk = YlxeEk 
the distance W[pk,o') > Ij^/k. 

To obtain a lower bound for W{pk,o') we use the dual formulation of 
the Kantorovich-Wasserstein distance Q and the function f{x) = d{x,Ek), 
which is in Lip^ Since d{x,Ek) = 0 on the support of pk we obtain 

W{pk,cr)> / d{x,Ek)da. 

Jx 

Vitali’s covering lemma ensures that for each k and for some e small enough, 
independent of k, there are at least 2ffEk pairwise disjoint balls of ra¬ 
dius ej^fk. Since the number of balls is twice the number of points in 
Ek, at least half the balls contain no point of Ek. We consider one such 
ball, B{xi,e/y/k). In the smaller ball B{xi,0.5e/Vk) we have d{x,Ek) > 
O.Be/Vk. Thus 


[ d{x,Ek)da [ d{x,Ek)da >'^-^a{B{xi,e/Vk)) 

Jx ,JB{xi,e/Vk) , Vk 

1 




1 

y/k 


Proof of Theorem To prove this we provide a (complex) transport plan 
between the probability measure hk{x) = ■^Kk{x, x) - bk stands for Bergman 
measure - and the empirical measure pk- We are going to prove that 

EW{pk,bk) = o(^^) . 
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Standard estimates for the Bergman kernel provide: 

hi. 


W{bt,u) = o[ 




Actually one can prove that the total variation (which dominates the Kantorovich- 
Wasserstein distance) satisfies: 


( 6 ) 


Kk{x,x) 


Nk 


— V 


< 


y/k 


This follows for instance by the expansion in powers of 1/A: of the Bergman 
kernel. In this context this is due to Tian, Gatlin and Zelditch, |Tia90 


Cat99,Zel98 


In the particular case of the spherical ensemble, the kernel is explicit and 
invariant under rotations, and the estimate is even better: the Bergman 
measnre is the equilibrium measure, i.e. bk = i'- 
Consider the transport plan 

^ Nk 

Pix,y) = — ^6xj{y){Kk{x,Xj),ej{x))da{x). 


i=i 


It has the correct marginals -- and respectively - and thns 

W{bk,yk)< // d{x,y)d\p\{x,y) 

J JXxX 

1 f 

<—V / d{x,Xj)\£j{x)\\Kk{x,Xj)\da{x). 

Now, letting /3' be the conjugate exponent of (3 (so that 1//3 + 1//3' = 1), we 
have 


(EIT)^ < 


< 


Ix^k 


'll yj^d{x, Xj)\£j{x)\\Kk{x, Xj)\da{x) j pk{xi,.. .,XNk 


(^j^d{x,Xj)\Kk{x,Xj)^ \(.j{x)\d\Kk{x,Xj)\d{x,Xj^ Pk{xi, 

Assnme for the moment that the following off-diagonal decay of the Bergman 
kernel holds: 

f C 

sup / d{x,y)\Kk{x,y)\da{x) < 


(7) 

yex Jx 
Then, by ([^, we obtain: 

P/P' r 1 Nk 


y/k 


f C \^ ^ f 1 f 

{EWf < j^\^j{x)\d\Kk{x,Xj)\d{x,Xj)pk{xl,... ,Xj,... ,XNk) 

(J \ 1^1 d' r \ r 

71 j Jx^k \^'^^^^^3)\dix,Xj)pk{xi,. . . ,X, . . . ,XNk)- 


,XNk)- 
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Finally, integrating first in Xj and applying again 0 we obtain 




Vk' 


"/k-' 


as desired. 

Estimate Q follows from the pointwise estimate for the Bergman kernel 

( 8 ) \Kk{x,y)\<CNke-^^‘^^^’y\ 


which holds when the line bundle is positive, see Ber03 


Indeed, consider the function h{s) = se strictly decreasing in +oo) 

For any y ^ X we bound the integral in Q as 

/• r+oo 

/ d{x,y)\Kk{x,y)\da{x) < / a {{x e X : h{d{x, y)) > s}) ds 
Jx Jo 

r+oo 1 

<Nk \h\s)\a{{x G X : d{x,y) < s}) ds < -j=, 

J{cVk)-^ vk 

where the last estimate follows from a{B{y, s)) < and ~ fe”. 

In the particular case of the spherical ensemble, the kernel is explicit and 
the decay is even faster: 


\Kkiz,w)\‘^ = [1- 


z — w\ 


(1 + |2:p)(l + |rcp) 


fc-i 


< Kk^ exp ( —Ck 


z — w 


' (1 + |7;p)(l + |tcP) 
= Kk^ exp [—Ck d{z, , 
where here d{z, w) coincides with the chordal metric. 


□ 


3. The determinantal setting 

Now we turn our attention to the almost sure convergence of the empirical 
measure. Using the fact that Lipschitz functionals of determinantal process 
concentrate the measure around the mean we prove the following result. 

Corollary 2. If yt is the empirical measure associated with the determi¬ 
nantal point process given by 0 with (3 = 2, and v denotes the equilibrium 
measure, then 

- //dimc(M) > 1 then IU(/ifc,i/) = 0{1/Vk) almost surely. 

- //dimc(M) = 1 then = 0{logk/\/k) almost surely. 

In particular, any realization of the spherieal ennsemble satisfies W{fiki = 
0{\ogk/y/k) almost surely. 

Let u be, as before, the normalized equilibrium measure. Let us define 
the functional / on the set of measures of the form a = ^xi by 

f{a) = nw(^^,uy 

As the Kantorovich-Wasserstein distance is controlled by the total variation, 
/ is a Lipschitz functional with Lipschitz norm one with respect to the total 
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variation distance. Here we use the following result of Pemantle and Peres 
(PPl4 , Theorem 3.5]. 

Theorem (Pemantle-Peres). Let Z he a determinantal point process of N 
points. Let f be a Lipschitz-1 functional defined in the set of finite counting 
measures (with respect to the total variation distance). Then 

F(f — E/ > a) < 3 exp ( —•—^-- 

J - ) - 16{a + 2N) 

Take now a = lOakNk/Vk, where ak = Vlog k for n = 1 and Ofc = 1 for 
n > 1. Then 

<F(^NkW{fik,i^)> NkEW{^rk,^i) + Wak^) 

lOOalN^/k 


,v) > 


VkJ 


< 3 exp — 


l6{l0akNk/Vk + 2Nk) 


< 


1 


ewi-^k^k/k) < p. 


Finally, a standard application of the Borel-Cantelli lemma shows that, 
with probability one, 


W{jik,v) < 


y/k 
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